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PREFACE.  This  technical  report  contains  a progress  summary  of 
recent  activities  In  IC  process  modeling  of  silicon  epitaxy  and 
oxidation.  The  report  builds  from  the  most  fundamental  process 
models  and  ends  with  details  of  work  In  progress.  Further 

details  appear  in  technical  reports  under  DARPA  contracts  DAA-B07- 
75-C-1344  «nd  DAA-B07-77-C-2684 . 

V 

ABSTRACT.  The  first  order  process  models  for  silicon  epitaxy  and 
oxidation  are  described.  Epitaxial  dopant  Inclusion,  autodoping 
and  transient  effects  are  discussed,  and  experimental  results  are 
presented.  Silicon  orientation,  surface  doping,  and  ambient 
effects  are  considered  for  slllcon-oxldatlon  rates.  . 
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1.  INTRODUCTION 


Epitaxy  and  oxidation  are  two  critical  steps  In  IC  technology. 
Process  models  date  back  to  the  mld-1960's.  Epitaxy  often  controls 
bipolar  circuit  performance  and  yields.  Oxide  growth  Is  essential 
In  all  IC  technologies,  and  segregation  frequently  dominates 
fabrication  limits  and  the  process  control  of  surface  effects.  The 
following  sections  will  review  the  basic  models  and  Introduce 
considerations  for  more  recent  experimental  work. 


SILICON  KPITAXY 


! 


2.1  Introduction 

I'lioro  is  an  extensive  literature  on  the  chemistry  itnd  reactor  dy- 
ii.imics  of  silicon  epitaxy  [ll  and  polyc  rys  tal  1 i ne  deposition  [2l. 

The  recent  advent  of  low-pressure  CVO  systems  will  undoubtedly  add 
substantially  to  this  body  of  knowledge  [3l. 

The  focus  of  this  discussion  is  on  process  modeling  for  silicon 
epitaxy  in  a horizontal  reactor  [ <1  1 . This  emphasis  is  based  on  the 
experience  at  Stanford  University  and,  in  large  part,  on  that  of 
technologists  using  commercial  equipment.  Figure  1 is  a schematic 
of  a typical  reactor.  Table  1 describes  the  various  source  gases 
and  their  functions,  and  Table  2 summarizes  the  key  features  of 
several  silicon  gaseous  sources.  Table  3 presents  a typical  re- 
actor sequence  used  during  epitaxial  silicon  growth.  The  follow- 
ing sections  will  elaborate  on  the  lundamental  models  that  describe 
the  kinetics  of  silicon  growth,  dopant  inclusion,  and  transient 
system  response. 

1. 


iig.  1.  Schematic  of  the  horizontal  epitaxial  reactor. 


Table  1 


GASES  USED  IN  EPITAXIAL  DEPOSITION 


i 


i 

[ 


Gas 

Function 

Comments 

N2 

main  flow 

Purges  out  ernloslve/polsonous  gases 
prior  to  opening  the  reactor  tube  to 
air 

»2 

main  flow 

Most  common  ambient  for  growth  or  epi- 
taxial layers 

SiCl, 

s 

SI  source 

Common  liquid  Si  source 
Vaporized  in  an  H2  bubbler 
Corrosive  vapor 

SiH, 

A 

Si  source 

Common  gaseous  Si  source 
Pyrophoric  gas 

HCL 

Si  etchant 

Most  common  Si  etchant  used  for  substrate 
preparat ion 
Corrosive  poison  gas 

PH3 

SI  dopant 

Most  common  phosphorous  source  for  doping 
epitaxial  silicon 
Flammable  poison  gas 

A8H3 

SI  dopant 

As  PH- 
3 

Sb(CH3)3 

SI  dopant 

A liquid  antimony  source  used  as  a vapor 
at  a concentration  of  a few  hundred 
ppm  in  H2 

Used  because  SbH3  is  unstable 
Poisonous  vapor 

®2»6 

Si  dopant 

As  PH3 

Table-  2 


SOURCES  OF  SILiaiN  IN  EPITAXY 


Typical  Conditions 

Source 

Temperature 

Rate 

Cosnents 

Reference 

rc) 

(p/mln) 

SIH, 

1000  to  10 SO 

0.2  to  1.0 

Pure  gaseous  source 

1,3, 4, 5. 8 

Pyrophoric  gas 
Low-temperature  deposition 

9.10,12 

Low  autodoping 

Moderate  growth  rates 

Surface  quality  sensitive  to  O2 

1 

13 

SlCl. 

1150  to  1200 

0.5  to  1.5 

Corrosive  liquid  source 

1,2,6,11 

Hlgh-teraperature  deposition 
Moderate  autodoping,  outd if fusion 
Moderate-high  growth  rates 
Most  common  source  for  linear 

13,14 

bipolar  integrated  circuits 
Easy  to  obtain  good  crystal 
quality  on  thick  layers 

15 

SlHCl. 

1150  to  1200 

1.0  to  10 

Corrosive  liquid  source 

1,7 

Hlgh-temperature  deposition 
Moderate  autodoping 
Very  high  growth  rates 
Most  common  source  of  poly-Si 
d lelect  r Ic  Isolation 

20 

Very  high  purity  epitaxial  layers 
used  in  high-voltage  dev/lces 

16,17 

SlHjClj 

1050  to  1100 

>1.0 

Caseous  source  at  7 psl 
Properties:  Intermediate  SiCl^p 

SIH4 

1,37 

SlBr^ 

1 

1 

1 

j 

— 

Rarely  used  source 

18 

Sl(CHj)^ 

1150 

0.4 

Rarely  used  source 

19 

Table  3 


TYPICAL  EPITAXIAL  GROWTH  CYCLE 


Step 

1 

Time 

Temperature 

G3  s 

Concentration 

Comment 

purge 

2' 

R.T. 

— 

Purge  out  O2 

H2  purge 

2' 

R.T. 

— 

Charge  to  H2  ambient 

Heat 

V 

1200°C 

— 

In  H2  ambient 

HCL  etch 

V 

1200°C 

1% 

Etch  0.26  u of  silicon 

H2  purge 

V 

to  1050°C 

— 

To  lower  temperature 
and  remove  HCL 

Growth 

00 

0 

1050°C 

0.05%  S1H4 
0.3  ppb  PH^ 

Growth  of  5 U,  1 fi-cm 
P-doped  Si  layer 

H2  purge 

r 

1050°C 

— 

Purge  reactants  prior 
to  cooling 

Cool 

4* 

R.T. 

— 

In  H2  ambient 

N2  purge 

2' 

R.T. 

— 

Before  opening  to  air 

2.2  Film  growth 

For  the  kinetics  of  film  growth,  the  silicon  source  gas  reacts  and 
decomposes  as  described  by  a typical  equation  for  silane, 

1000“C 

SIH,  SI  + 2H„  (1) 

4 2 

The  silicon  growth  rate  is  proportional  to  the  partial  pressure  of 
silane , 


V = Kj(SiH^)j 


(2) 


where  V Is  the  growth  rate,  Kx  is  the  surface-reaction  rate  con- 
stant, and  the  subscript  denotes  the  surface  quantities.  For 
halide-bearing  source  gases,  there  are  reactive  conditions  whereby 
silicon  etching  occurs;  however,  for  model  simplicity,  these 
source  gases  will  not  be  considered  because  of  the  complexity  of 
the  multispecies  reaction. 

Equation  (2)  describes  growth  only  under  the  conditions 
of  surface-reaction  control;  that  is,  surface  reaction  is  the 
limiting  step.  Figure  2 is  a typical  Arrhenius  plot  for  silane 
growth  vs  Inverse  temperature  with  gas  partial  pressure  as  a 
parameter  L4] . The  silane  concentrations  are  those  of  the 
turbulent  gas  stream. 


TEMPERATURE  (°C) 

1 100  1000  900  800 


Fig.  2.  Arrhenius  plot  of  SiH^  growth  rate. 
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Figure  3 Is  a cross  section  of  the  reactor.  A boundary 
layer  [5]  (often  labeled  a "stagnant"  layer)  exists  between  the 
turbulent  gas  and  reactive  wafer  interface  T;  also  shown  are 
plots  of  reactant  gas  species  under  two  limiting  cases. 
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Fig.  3.  Cross  section  of  horizontal  reactor  with  thermal  and 
chemical  profiles  during  growth. 


Epitaxial  growth  proceeds  by  the  following  steps: 

(1)  mass  transfer  of  the  reactant  molecules  (such  as  SiH^) 
by  diffusion  from  the  turbulent-layer  reservoir  across 
the  boundary  layer  to  the  silicon  surface. 

(2)  adsorption  of  reactant  atoms  on  the  surface 

(3)  one  or  more  chemical  reactions  at  the  surface 

(4)  desorption  of  product  molecules  (such  as  H^) 

(5)  mass  transfer  of  the  product  molecules  by  diffusion 
through  the  boundary  layer,  back  to  the  turbulent 
layer 

(6)  lattice  arrangement  of  tlie  adsorbed  silicon  atoms 
(may  occur  as  part  of  3) 

Concentration  gradients  across  the  boundary  layer  are  such  that 
there  Is  a diffusion  flux  of  product  molecules  away  from  the 
surface.  The  turbulent  layer  is  depicted  as  having  no  concentra- 
tion or  temperature  gradient. 


The  overall  deposition  rate  is  determined  by  the  slowest 
of  processes  (1  to  b above).  The  expected  temperature  depend- 
ence of  the  growth  rate  varies  markedly  for  different  controlling 
mechanisms.  If  diffusion  of  either  reactant  or  product  across 
the  stagnant  layer  (so-called  mass-transport  control)  is  the 
slowest  part  of  the  ractlon,  the  growth  rate  would  not  depend 
on  the  temperature  to  first  order.  On  the  other  hand,  if  a 
surface  chemical  reaction  Is  the  slowest  process  (surface-reaction 
rate  control),  one  would  expect  the  deposition  rate  to  have  the 
same  temperature  dependence  as  the  chemical  reaction,  which  is 
an  exponential  function  of  Inverse  temperature.  Furthermore, 
because  the  various  silicon  sources  undergo  different  chemical 
reactions  with  different  activation  energies,  one  would  expect 
them  to  have  diverse  growth-rate  temperature  dependences  under 
surface-reaction  rate  control. 


The  cu  'wn  in  Figure  2 displays  a low-temperature 

reaction  r which  surface-reaction  rates  dominate  the 

growth  r .igh-temperature  nonactivated  reaction  realm 

in  whlc'  ; ansport  dominates  the  grovTih  rate.  One  striking 

feature  is  tnat  similar  plots  for  SiCl^,  SIH2CI2,  and 

SiHCli  have  the  same  activation  energy;  (roughly  3/  kcal/mole  or 
1.6  eV/molecule) . This  common  activation  energy  in  the  temperature- 
dependent  region  implies  that  the  dominant  surface  "reaction"  is 
step  6;  (the  lattice  arrangement  of  the  silicon  atoms).  Indeed, 
the  activation  energy  for  surface  diffusion  of  adsorbed  silicon 
atoms  on  silicon  has  been  measured  to  be  36  1 6 kcal/mole  [6]. 


At  higher  temperatures,  growth  reaction  is  nonactivated, 
suggesting  that  the  slowest  process  is  mass  transport  of  either 
products  or  reactants.  The  rate-limiting  step  is  probably  react- 
ant transfer.  This  hypothesis  is  supported  by  two  considerations. 
First,  in  the  boundary  layer  under  reactant  mass-transport  limited 
conditions,  the  growth  rate  is  proportional  to  the  diffusion  co- 
efficient of  the  reactant  species  in  hydrogen.  This  coefficient, 
in  turn,  is  proportional  to  the  inverse  square  root  of  the  molecular 
weight  of  the  reactant  molecule.  Second,  the  large  temperature 
gradient  in  the  boundary  layer  tends  to  retard  the  gaseous  diffusion 
of  reactants  toward  the  susceptor  and  to  increase  the  diffusion  of 
products  away  from  the  susceptor. 

In  summary,  boundary-layer  theory  and  empirical  observations 
produce  the  following  qualitative  results  for  the  epitaxial  growth 
rate  of  common  silicon  sources. 

(1)  The  growth  rate  indicates  two  distinct  regions,  a low- 
temperature  region  in  which  the  growth  rate  fits  an 
Arrhenius  plot  and  a high-temperature  regi'. .1  in  which 
the  growth  rate  does  not  depend  on  temperature. 


(2)  In  the  low- temperature  region: 

(a)  Activation  energy  is  approximately  37  kcal/mole, 
a number  consistent  with  the  activation  energy 
for  surface  diffusion  of  adsorbed  silicon  atoms. 

(b)  Activation  energy  is  independent  of  source  type. 

(c)  For  low  source  concentrations,  the  growth  rate  is 
directly  proportional  to  the  input  silicon  source 
concentration. 

(d)  The  linear  growth  rate/source  concentration  rela- 
tionship applies  at  high  growth  rates  only  for 
SiH^.  The  cholorlde-containing  sources  do  not  obey 
this  relationship  because  a reverse  reaction,  HCl 
etching,  occurs  [7]. 

(3)  In  the  high-temperature  region: 

(a)  The  temperature  Insensitive  growth  rate  is  con- 
trolled by  mass  transport  of  the  silicon  source 
from  the  turbulent  layer  to  the  silicon  interface 
by  diffusion. 

(b)  For  low  source  concentrations,  the  growth  rate  is 
directly  proportional  to  the  input  silicon  source 
concentration. 

(c)  The  growth  rates  for  chloride-containing  sources 
are  not  linear  functions  of  source  concentration 
at  high  concentrations  because  of  HCl  etching. 

The  changeover  between  the  surface-controlled  and  mass- 
transport  controlled  realms  of  SiH,  deposition  can  be  described  by 
the  balancing  of  molecular  fluxes  because  no  reverse  reactions 
complicate  the  SIH^  problem.  The  growth  rate  V is  directly  pro- 
portional to  the  SiH4  concentration  at  the  silicon  surface,  as 
given  by  Eq.  (2).  The  surface-reaction  rate  constant  K has  the 
characteristic  activation  energy  of  37  kcal/mole  and  the  units  of 
the  growth  rate  (cm/sec)  because  [SlH^Jj  is  a unitless  ratio  of 
gas  flows.  The  flux  of  SiH^  atoms  across  the  boundary  layer  F^ 
is  proportional  to  the  SlH^^concentratlon  gradient  across  the 
layer,  which  is  the  difference  between  SIH^  concentration  in  the 
turbulent  layer  [SiH^]^  and  that  at  the  Interface  [SiH^]^, 

Fj  - h(SlH^)|[SiH^].^  - [SiH^]j|  (3) 

where  h(SlH^)  is  the  mass-transport  coefficient  for  SiH^  and  is  a 
temperature- insensitive  constant  that  relates  molecular  flux  to 
the  concentration  gradient;  it  is  proportional  to  the  SIH^  diffusion 
coefficient  in  H2  and  is  Inversely  proportional  to  the  boundary- 
layer  thickness.  The  flux  of  silicon  atoms  being  Incorporated  into 
the  lattice  is 


(4) 


^’2  = "Sl" 


22  -3 

in  which  N„.  is  the  atomic  deii'^ity  of  silicon,  N . = 5.0  x 10  cm 
A simultaneous  solution  of  Eqs.  (2)  through  (4)  oDtalns  the  inter- 
*^ace  silane  concentration, 


Si 


h(SlH, ) 
4 


There  are  two  limiting  cases; 


(5) 


Case  1:  Surface-reaction  control. 

At  low  temperatures,  K is  small  enough  so  that  [SiH, ] [SiH, ] 
(See  Figure  3).  In  this  case,  from  Eq.  (2), 

V = Kj[SiH^]^  (6) 


Case  2:  Mass-transport  control. 

At  high  temperatures,  K becomes  large,  and  [SiH,]  = 0;  then 
from  Eqs.  (2)  and  (3), 


h(SlH, ) 

''  ■ <') 

The  rate  constant  for  mass-transport  control  is  K^(cm/sec). 

Although  Eqs.  (6)  and  (7)  predict  a first-order  growth-rate 
dependence  on  silane  concentration,  the  surface-reaction  rate 
constant  is  exponentially  dependent  on  temperature  and  the  mass- 
transport  rate  constant  has  only  a linear  dependence  on  tem- 
perature (essentially  no  dependence  over  the  narrow  1000°  to 
1100°C  range  commonly  used  for  SiH^  epitaxy) . 

Experimental  SIH^  growth  rate  is  plotted  vs  temperature  in 
Figure  2.  The  surface-reaction  region  has  an  activation  energy  of 
AE  * 39  kcal/mole.  Transfer  to  the  mass-transport  region  is 
nearly  complete  at  1000°C  for  all  growth  rates  from  0.1  to  0.5  u/mln. 
From  1000°C,  the  maximum  variation  observed  is  only  10%.  Figure  4 
is  a graph  of  growth  rate  vs  silane  concentration  for  T “ 1050°C. 

The  slope  of  t^e  line  is  the  mass-transport  rate  constant  value 
' 7.5  X 10  cm/sec  = 450  y/min.  These  values  for  activation 
energy  [6]  and  the  rate  constant  [5]  are  consistent  with  the 
literature. 
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[SiH4]^  = FLOW  (SiH^l/FLOW  (Hg) 


Fig.  4.  SiH^  growth  rate  vs  SIH^  concentration. 

2.3  Doping  of  epitaxial  layers 

2.3.1  Injected  species.  Epitaxial  layers  are  most  commonly 
doped  by  the  Incorporation  of  a small  amount  of  dopant  hydrides 
(PH-,  AsH-,  B^H^) . Phosphine,  arsine,  and  dlborane  are  gaseous 
sources  that  can  be  diluted  to  the  10“^  level  required  for  lightly 
doped  epitaxial  layers.  The  doping  of  SIH,  epitaxial  layers  with 
dopant  hydride  sources  has  been  described  theoretically  [8,9]. 

The  dopant  Is  transported  across  the  boundary  layer  by  diffusion 
as  described  In  the  previous  section  for  SIH4.  Tne  analysis  Is 
complicated,  however,  by  the  fact  that  the  Incorporated  dopant  Is 
nonunlformly  distributed  near  the  surface  of  the  growing  layer. 

The  relevant  dopant  concentrations  In  the  reactor  and  In  the  wafer 
are  presented  In  Figure  5 for  PH3.  The  turbulent  layer  (Input) 
dopant  molar  concentration  [PH3)'j'  Is  the  flow  ratio  PH3/H2.  The 
concentration  of  phosphine  near  the  surface  of  the  growing  epitaxial 
layer  Is  [PH3]j.  The  concentration  of  Incorporated  phosphorous 
atoms  near  the  surface  of  the  layer  [P]j  Is  higher  than  the 
phosphorous  concentration  deeper  In  the  bulk  material  [Pig  [f.9j. 

All  of  the  Incorporated  phosphorus  Is  considered  to  be  Ionized. 

The  solid-state  phosphorous  concentrations  [Plj  and  [Pjg  have  units 
of  (atoms  P/atoms  SI)  and  must  be  multiplied  by  the  atomic  density 
of  silicon  “ 5.0  X 10^2  atoms/cm^)  to  obtain  the  conventional 

doping  concentration. 
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Fig.  5.  Dopant  concent  rat  ions  in  the  reactor  anti  in  the  wafer. 

(a)  Dopant  concentrations  near  an  epitaxial  Interface 

(b)  Band  bending  near  a growing  silicon  surface 

(c)  Dopant  and  free-electron  concentrations  near  a 
growing  silicon  surface 


The  goal  of  this  analysis  is  to  predict  the  solid-state 
doping  [P]g  from  the  turbulent-layer  phosphine  concentration  [PHgJ-j- 
as  a function  of  temperature  and  growth  rate.  The  effective 
segregation  coefficient  relates  these  two  variables  as 


[P], 


[PH 


3^T 


EFF 


(8) 


For  doping  concentrations  lower  than  the  concentration  of  free 
electrons  in  silicon  at  epitaxial  temperatures  ( 10^"  cm~^), 
the  segregation  coefficient  is  a function  of  only  tempevature  and 
growth  rate,  but  not  doping  level  [8];  that  is,  for  a set  temperature 
and  growth  rate,  solid-state  doping  is  a linear  function  of  gas- 
phase  dopant  concentration.  With  respect  to  the  dopant  concentra- 
tions defined  in  Figure  5a,  Eq.  (8)  can  be  expressed  as 


f^^B  1 

[p]b 

[Pll 

[PHslj 

[PH  ] Y 

E r 

[P]j  • 

[PH3]j  • 

(a) 

(b) 

(c) 

The  overall  partition  coefficient  Ygpp  can  now  be  evaluated  term 
by  term.  Ratio  (a)  on  the  RHS  of  Eq.  (9)  relates  the  surface  and 
bulk  solid-phase  dopant  levels  that  may  be  ditferent  because  of 
band  bending  at  the  growing  silicon  surface.  Ratio  (b)  relates 
the  actual  solid  and  gas  doping  concentrations  at  the  Interface. 
Ratio  (c)  is  essentially  a measure  of  the  boundary  layer  doping 
gas  gradient  which  is  the  driving  force  that  provides  the  required 
flux  of  dopant  to  the  surface. 


Following  the  notation  in  Figure  5 and  the  derivations  of 
Rodgers  [4]  and  Bloem  [8],  the  expression  for  an  effective  segre- 
gation coefficient  is 


nil 


J; 

^EFF 


exp  — + exp-^  1^1  - exp 


h(SiH^) 

where  (})„  is  the  surface  potential  pinned  at  one-third  the  gap 
above  tne  valence  band  [8],  <|)g  la  the  bulk  Fermi  potential,  and 
is  a critical  growth  rate  at  which  trapping  starts  to  occur. 

The  parameter  h(PH.)  is  the  phosphine  transport  coefficient  similar 
to  h(SlH^),  and  the  other  parameters  have  their  standard  meanings. 
Figure  5 aefines  the  equilibrium  segregation  term 


i ^ 

Y " [PH3]j 


. -1 


[siH^J-r 


(10) 


(11) 


For  low  growth  rates,  the  effective  segregation  coefficient  Is 
equal  to  the  equilibrium  segregation  coefficient  Y»  modified  to 
account  for  trapping.  Reactant  Incorporation  Is  surface-reaction 
limited.  For  high  growth  rates,  the  effective  segregation  coef- 
ficient Is  proportional  to  the  Input  silane  concentration  [SIH^]^. 
Four  limiting  cases  of  Eq . (10)  exist  as  a result  of  combining 
the  two  limiting  conditions  in  the  crystal  (bulk-surface  equilibrium 
and  complete  trapping)  with  the  two  limits  for  dopant  transport  to 
the  crystal  (surface  reaction  and  mass  transport).  These  cases  are 
depicted  schematically  In  Figure  6 for  an  arbitrary  donor  hydride 
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Fig.  6.  Limiting-case  dopant-concentration  profiles  for  silicon 
epitaxial  layers. 


Case  1. 


Surface-reaction  limited  dopant  transport  and  bulk-surface 
equilibrium  conditions  would  occur  during  growth  at  moderately 
high  temperatures,  low  growth  rates,  and  with  a fast  diffusing 
dopant  (1150°C,  0.1  y/mln,  PH3)  [9].  In  this  case,  Eq.  (10)  can 
be  approximated  by 
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The  effective  segregation  coefficient  is  the  equilibrium  segregation 
coefficient  modified  by  band-bending  effects. 

Case  2. 

Surface-reaction  limited  dopant  transport  and  surface  dopant 
trapping  conditions  would  occur  during  growth  at  low  temperatures, 
moderate  growth  rates,  and  with  a slow  diffusing  dopant  (1000°C, 

0.5  y/mln,  AsH_).  The  effective  segregation  coefficient  is  then 
the  equilibrium  segregation  coefficient. 
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Case  3. 

Mass-transport  limited  dopant  delivery  and  bulk-surface 
equilibrium  conditions  would  occur  during  growth  at  very  high 
temperatures,  moderately  high  growth  rates,  and  with  a fast  diffusing 
dopant  (1225°C,  2 y/mln,  PH3).  The  effective  segregation  coefficient 
is  determined  by  mass  transport  of  dopant  across  the  boundary 
layer.  Equation  (10)  can  be  approximated  as 

_ _1 _ 1 

[PHsIt  “ Yeff  ’ [h(SiH^)/h(PH^)]  ^ 

Case  A . 

Mass-transport  limited  dopant  delivery  and  surface  dopant 
trapping  conditions  would  occur  during  growth  at  relatively  low 
temperatures,  high  growth  rates,  and  with  a slow  diffusion  dopant 
(1050°C,  4 y/mln,  ASH3).  The  effective  segregation  coefficient  is 
given  in  Eq . (14) . 
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Although  all  ol  th''  growth  conditions  depicted  in  Fig.  6 are 
possible,  practical  con.s  iderat  ions  limit  Sill4  epitaxial  growth  to 
the  1000°  to  1050°C  range  and  a 0.2  to  1.0  p/min  growth  rate.  At 
higher  temperatures,  the  SIH^  process  has  the  following  problems. 

(1)  At  high  growth  rates,  S1H4  tends  to  decompose  in  the 
boundary  layer,  forming  a "cloud"  in  the  reactor 
that  quickly  coats  the  exit  end  of  the  reactor 

tube . 

(2)  Silicon  deposition  on  the  reactor  tube  Increases 
for  air-cooled  tubes. 

(3)  Autodoping  and  outdif fusion  begin  to  become  problems. 

At  temperatures  below  i000°C,  crystal  quality  suffers;  there 
is  a maximum  epitaxial  growth  rate  for  evey  temperature  below 
which  polycrystalline  deposition  occurs.  Growth  rates  below  0.2 
p/min  are  impractical  for  common  3 to  20  p layers  because  long 
deposition  times  are  not  only  wasteful,  but  lead  to  excessive  out- 
diffusion.  Typical  silane-reaction  conditions  are  intermediate, 
therefore,  to  the  conditions  depicted  in  Figure  6a  and  6b. 

Some  surface-dopant  trapping  is  present;  an  increase  in  growth 
rate  leads  to  an  increase  in  solid  doping  level,  all  other  variables 
held  constant. 

The  effective  segregation  coelf  icient  has  been  measured  ex- 
perimentally over  a wide  range  ol  conditions.  Bulk  doping  concen- 
tration vs  monotomlc  gas-phase  dopant  concentration  is  presented 
in  Figure  7 lor  a combination  of  different  dopant  s(jurces,  sili- 
con sources,  growth  rates,  and  temperatures.  The  doping  concen- 
tration is  calculated  as  the  effective  flow  ratio  of  dopant  gas 
to  main  hydrogen  flow.  With  reference  to  the  reactor  in  Figure 
1,  this  ratio  is 


I f (dopant) . f (mix . f ( Inject) 

3 f(dopant‘  i £(lst  H^)  f (mix)  + f(2nd  H2)  f(main) 

(15) 

where  f is  the  flow  in  liters  per  minute.  The  concentration  in 
the  right-hand  scale  is  the  normal  concentration  per  cubic  centi- 
meter divided  by  All  the  dopants  have  higher  concentrations 

in  the  bulk  than  in  the  gas  phase  by  a factor  of  10  to  200  and 
have  approximately  the  same  effective  partition  coefficient 

(«=0.01)  except  for  AsH.,/SlH,  at  1050°C  and  (CH.J , Sb/SIH,  . 

> 4 3 3 4 

Tlie  fact  that  the  effective  partition  coefficients  are 

independent  of  growth  rate,  temperature,  and  even  silicon  source 

for  PH  and  B H suggests  a surface-reaction  limited  process  with 
> 2 o 
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Fig.  7.  Solid-state  doping  vs  gas-phase  dopant  concentration 

for  various  dopant  sources,  silicon  sources,  and  temperatures. 


bulk-surface  equilibrium  (Figure  6(a)).  Trapping  is  not  present 
to  any  appreciable  degree  for  PH^/SiH^  because  of  the  low  0.25 
p/min  growth  rate.  For  PH3/S1C14,  821(5/31114,  the  high  growth 
temperatures  (1075°  to  1180°C)  establish  surface-bulk  equilibrium. 
There  is  some  trapping  in  the  1050°C  ASH3/SIH  system  relative  to 
the  1180°C  AsH^/SiCl^  system,  as  evidenced  by  the  higher  solid 
doping  concentration  for  a given  input  arsine  concentration. 

The  trapping  phenomenon  in  ASH3/SIH4  epitaxial  layers  has 
been  considered  as  a function  of  temperature  and  growth  rate. 
According  to  Eq.  (10),  doping  should  decrease  as  temperature  (V  ) 
Increases  in  the  presence  of  trapping.  This  trend  is  verified*" 
for  three  growth  rates  in  Figure  8.  All  the  curves  begin  to 
flatten  out  at  1025°C,  indicating  that  trapping  is  nearly  complete, 
(Aslg  « (As)^. 
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Fig.  S.  Resistivity  of  arsine-doped  silane-grown  epitaxial 
layers  as  a function  of  temperature  with  growth  rate  and 
arsine  concentrations  as  parameters. 


2.3.2  Autodoping.  One  of  the  limiting  and  most  important 
considerations  in  epitaxial  deposition  is  the  movement  of  substrate 
impurities  into  epitaxial  layers.  Two  components  are  considered  — 
outdiffusion  (mt)vement  caust J by  solid-state  thermal  diffusion  in 
silicon)  and  autodoping  (extia  movement  of  substrate  doping  into 
an  epitaxial  layer,  which  cannot  be  explained  by  diffusion  effects). 
The  following  dopant  fluxes  from  a wafer  into  an  epitaxial  layer 
must  bi‘  considered: 

(1)  dopant  "evaporated"  from  the  back  of  the  wafer,  which 
mixes  into  tlie  turbulent  layer 

(2)  dopant  "evaporated"  from  the  front  of  the  wafer  into 
the  boundary  layer 

(3)  dopant  that  outdiffuses  from  the  substrate  ;is  a result 
of  therm.ll  diffusion  alone. 

For  standard  biirled-col lector  bipolar  processing,  a moderately 
doped  (10^^  to  10*^  cm"  ^)  epitaxi.il  layer  is  grown  over  localized 


n"*^  buried  layers.  'Hie  backside  of  tlie  wafi.r  is  usually  masked 
during  the  buried-layer  diffusion  to  eliminate  this  autodoping 
source.  Tlie  wafer  is  then  subjected  to  very  high  temperature- 
time cycles  to  produce  junction  isolation;  in  this  case,  the 
sealed-wafer  backside  does  not  provide  an  autodoping  source,  and 
frontside  autodoping  is  dominated  by  outdiffusion  from  the  high- 
temperature  isolation  step.  Consequently,  for  standard  buried- 
collector  processing,  outdiffusion  can  explain  dopant  migration 
♦"rom  the  burled  layer  into  the  epitaxial  layer  [I0|. 

For  low-temperature  silane  epitaxial  growth,  especially 
when  high-temperature  processing  steps  do  not  follow,  there  is  a 
definite  movement  of  dopant  into  the  growing  layer  from  the  substrate 
that  cannot  be  explained  by  outdiffusion.  Figure  9 is  the  profile 
of  an  undoped  epitaxial  layer  grown  with  silane  (1100°C,  0.4  y/min) 
on  an  arsenic-doped  substrate.*  For  the  10  min  period  at  high 
temperature,  the  diffusion  length  for  arsenic  is  /DE  = 0.041  y. 

The  profile  solution  for  outdiffusion  [10]  is 

N(epi)  = ^ N(substrate)  erfc  / — - — \ (16) 

^ '2  /Dt' 

in  which  x is  the  distance  into  the  epitaxial  layer  measured  from 
the  Interface.  The  calculated  outdiffusion  curve  is  graphed  in 
Figure  9. 


DISTANCE  FROM  EPITAXIAL  INTERFACE  (/i) 

Fig,  9.  Autodoping  profile  for  an  intrinsic  silane-grown  epitaxial 
layer  on  an  arsenic-doped  substrate. 


*A11  profile  measurements  were  ntade  using  the  spreading-resistance 
probe  technique. 


In  this  case,  the  dop.inl  movement  (lure  ealled  ",i  iit  tuJop  i n^")  i s 
dominant  over  the  "out d i 1 1 us  ion"  term.  Ueeause,  for  this  experi- 
ment, the  arsenic-doped  "wafer"  was  aitually  an  a rsen i t -doped 
epitaxial  layer  deposited  on  a li^litly  boron-doped  water,  the 
sole  component  of  autodoping  contributing  to  the  measured  profile 
was  evaporation  from  the  front  of  t'ue  wafer  into  ttie  growing  layer. 

A well-known  feature  of  the  frontside  evaporation  solution  is 
that  the  evaporation  rate  is  proportional  to  the  surface  concen- 
tration [11].  As  the  surface  concentration,  resulting  from  auto- 
doping decreases,  the  evaporation  rate  decreases  which,  in  turn,  i 

reduces  concentration  in  the  growing  layer;  the  doping  profile 
is  a simple  decaying  exponential  extrapolating  back  to  an  effective 
surface  concentration  N^*  which  is  determined  empirically  for  a 
given  dopant,  silicon  source,  and  set  of  growth  conditions.  For 
arsenic-silane  in  Figure  9,  the  profile  is  close  to  exponential, 
dropping  a factor  of  10  every  2.67  p. 

This  behavior  has  also  been  measured  in  layers  grown  ai 
higher  temperatures  with  SiCl4.  Figure  10  shows  the  measured 
autodoping  tail  of  a boron-doped  substrate  into  an  undoped  epitaxial 
layer  grown  with  SiCl^  at  1180°C  and  1.0  y/min.  Here,  the  calculated 
outdif fusion  is  comparable  to  the  autodoping.  Outdif fusion  is  much 
more  prominent  than  in  the  previous  case  because  of  the  higher 
growth  temperature  and  the  higher  diffusivity  of  boron  relative  to 
arsenic.  Again,  an  exponential  fit  is  within  experimental  error. 

The  curve  decays  a factor  of  10  every  0.86  y.  .If  a substantial 
subsequent  heating  step  such  as  diffusion  isolation  had  been  done 
to  the  substrates  used  in  the  boron  SiCl^  experiments,  the  outdiff- 
usion  tail  would  have  become  dominant  over  autodoping. 


Fig.  10.  Autodoping  profile  fo 
for  an  intrinsic  silane-grown 
epitaxial  layer  on  a boron- 
doped  substrate. 


One  Immediate  conclusion  that  can  be  drawn  from  a comparison 
of  Figures  9 and  10  is  that  the  dopant  species  Is  a very  Important 
consideration.  One  would  expect  a priori  the  slow-diffusing 
dopant,  As,  grown  at  the  lower  temperature  to  have  much  less  auto- 
doplng  than  the  fast-dlf fusing  dopant,  B,  grown  at  the  higher 
temperature.  Because  the  opposite  Is  true,  arsenic  must  be  a 
poorer  autodoping  element  than  boron.  One  significant  factor  Is 
vapor  pressure.  Arsenic  "evaporates"  more  easily  than  boron. 

The  exact  nature  of  the  autodoping  differences  between  the 
elements  have  not  been  reported,  but  empirical  data  to  be  presented 
here  demonstrate  that  the  worst- to-best  autodoping  order  of  common 
silicon  dopants  is  As  >■>  B > P > Sb. 

The  second  source  of  dopant  for  autodoping  Is  evaporation 
from  the  backside  of  the  wafer.  The  amount  of  autodoping  Is  pro- 
portional to  the  evaporation  rate.  This  rate  and  the  autodoping 
profile  depend  to  first  order  on  wafer  dopant  concentration  N^. 

The  net  effect  Is  that  the  wafer  backside  supplies  dopant  to 
the  epitaxial  layer  at  a rate  Initially  set  by  Its  concentration, 
evaporation  rate,  dlffuslvlty,  and  temperature.  The  amount  of 
dopant  supplied  then  decreases  slowly  In  time  because  of  dopant- 
depletion  effects  on  the  wafer  backside.  The  time  variation  Is 
usually  weak  so  that  the  backside  appears  as  a relatively  constant 
autodoping  source.  The  solution  to  this  boundary  condition  Is 
simple.  For  a given  set  of  growth  conditions,  the  autodoping 
profile  Is  constant  and  dependent  only  on  Initial  wafer -backside 
doping.  In  a study  [12]  of  Intrinsic  epitaxial  layers  (SIH^, 

1050°C,  0.35  y/mln)  grown  on  arsenic-doped  wafers,  the  following 
empirical  relation  was  derived: 

N(constant  autodoping  level)  = (2.5  x 10^)  N(backslde)  (17) 

The  constant  autodoping  levels  were  measured  at  4 to  7u  from  the 
original  wafer  surface,  after  the  point  at  which  the  frontside 
evaporation  and  outdlf fusion  transients  had  died  out. 

With  a general  description  of  the  three  components  of  dopant 
movement  Into  an  epitaxial  layer  (outdlf fusion  and  frontside  and 
backside  autodoping),  a conceptual  picture  of  autodoping  can  be 
drawn.  Figure  11  shows  that  outdlffuslon  dominates  nearest  the 
Interface,  frontside  evaporation  effects  become  apparent  as  soon 
as  the  outdlffuslon  tall  dies  out,  and  backside  evaporation  dominates 
the  autodoping  profile  farthest  from  the  Interface  because  the 
backside  autodoping  source  remains  relatively  constant  and  the 
frontside  source  Is  continuously  burled  under  growing  silicon. 

Figure  11  Is  highly  Idealized;  for  example,  the  analysis  does  not 
consider  lateral  autodoping  variations  that  Include  systematic 
worse  autodoping  at  the  wafer  edges  as  a result  of  the  backside 
source.  Any  of  the  three  regions  In  Figure  11  may  not  be  present 


under  certain  conditions,  (large  outdif fusion  tail  may  go  completely 
to  the  surface  of  a thin  epitaxial  layer). 


DISTANCE  FROM  EPITAXIAL  INTERFACE  (/i) 

Fig, 11.  Conceptual  diagram  of  autodoping  and  outdiffuslon  profiles. 
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2.3.3  Time-dependent  effects.  The  above  subsections  have 
defined  injection  and  autodoping  effects  primarily  for  the  steady- 
state  condition.  Recent  work  has  demonstrated  practical  means  for 
obtaining  a predefined  doping  profile  by  utilizing  transient  doping 
effects.  A two-step  process  [13]  using  silane  has  achieved  abrupt 
epitaxial  Interfaces  by  allowing  the  boundary  layer  vapor  pressure 
Pg  to  decay  with  time  after  the  growth  of  a first  layer.  The  concen- 
tration in  the  first  layer  is  reduced  from  the  bulk  value  Ng  to  a 
value  set  by  the  ratio  of  loss  coefficients  Cgj  from  bulk  to  the 
first  boundary  layer,  and  from  the  second  boundary  layer  to 
the  epitaxial  layer, 


N * • — 
B K 


BI 


(18) 


Following  this  decrease  in  surface  concentration,  the  vapor  pressure 
decays  to  a second  equilibrium  value  with  an  exponential  time 
behavior.  Subsequent  epitaxial  growth  after  the  decay  of  Pg  results 
in  an  abrupt  doping  profile. 

A transfer-function  characterization  method  has  been  reported 
that  predicts  profiles  for  time-varying  injected  dopant  species 
[ill.  Figure  12a  is  a block  diagram  of  the  epitaxial  reactor  as  a 
system,  the  input  is  the  gas  flow  as  a function  of  time  f(t),  and 
the  output  is  the  doping  density  profile  in  the  epitaxial  layer 
N(x).  All  four  blocks  Interact,  With  respect  to  the  doping 
mechanism  and  despite  such  complexity,  the  reactor  can  be  considered 
analogous  to  a linear  system  within  limitations.  The  implication 
of  this  is  that  it  is  possible  to  obtain  a transfer  function  of 
the  reactor  that  relates  the  dopant  gas-flow  input  f(t)  to  the 
doping  profile  output  N(x)  as  shown  in  Figure  12b.  After  the 
transfer  function  has  been  determined,  it  is  possible  to  calculate 
the  dopant  gas  flow  as  a function  of  time  required  to  achieve  a 
desired  dopant  profile  in  an  epitaxial  film. 
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Fig,  12.  Block  diagram  of  reactor  system. 

(a)  Cross  section  of  the  reactor  chamber 

(b)  The  reactor  as  a system 

(c)  One  transfer  function  relating  the  dopant  gas-flow 
Input  f(t)  to  dopant  concentration  In  the  epitaxial 
film  N(x) 


Tniii,'-;  i(Mi  t I'l^sponse  oI  t luj  tlopani  incorporation 
mechanism  . A linear  system  can  lie  cliaracle  ri  zed 
l)V  ot)tainin^;  its  transient  response.  Alter  tliis  response  is 
known,  it  is  .t  relatively  tnisy  task  to  <iet(;rmine  the  transler 
1 line  t ion  . 

Although  there  are  definite  limits  to  this  approach,  it  has 
been  used  successfully  in  several  chemical  systems.  For  example, 
diffusion-rate  processes  in  catalytic  reactors  have  been  characterized 
by  "admittance  functions"  analogous  to  ac  electrical  parameters 
[15].  The  time-varying  output-gas  concentrations  were  related  to 
the  time  variations  of  the  input  gases.  Further  work  has  considered 
diffusion  through  a boundary  layer,  adsorption  on  a surface,  and 
diffusion  into  a solid  [16].  Several  of  the  same  kinetic  processes 
are  present  in  our  study  of  dopants  in  epitaxial  layers;  however, 
additional  chemical  processes  are  involved  because  the  output  is 
not  a time-varying  gas  concentration,  but  a variation  of  the  solid 
dopant  concentration  as  a function  of  position  in  the  deposited 
epitaxial  layer. 

In  addition  to  providing  a method  for  fabricating  a desired 
dopant  profile,  this  investigation  will  produce  basic  information 
concerning  the  various  mechanisms  involved  in  the  dopant-inclusion 
processes.  As  indicated  by  Kobayashi  and  Kobayashi  [17],  transient 
studies  may  reveal  the  importance  of  various  mechanisms  that  may 
not  be  apparent  in  steady-state  studies.  In  particular,  consider- 
ation of  limiting  cases  may  allow  separation  of  the  several  mechanisms 
Involved;  dopant  t;as  1 low  without  film  growth,  for  example,  would 
Involve  only  a selected  number  of  mechanisms  involved  in  the  total 
epitaxial  doping  process. 

The  reactor  was  first  optimized  for  a nominal  deposition  of 
approximately  0.6  p/mln  and  uniform  doping  during  the  entire 
deposition  process.  Silane  was  the  source  of  Si,  and  deposition 
was  accomplished  at  1070°C.  Arsine  was  selected  as  the  dopant 
gas,  and  flow  settings  were  determined  for  the  typical  partial 
pressure  of  approximately  2 x 10~  and  6 x 10“^  atm,  corresponding 
to  dopant  concentrations  of  roughly  1 x 10^^  and  3 x 10  cm“^  in 
the  epitaxial  layer.  The  substrates  were  <100>-oriented. silicon 
wafers  with  phosphorous  dopant  concentrations  in  the  10^  cm“^ 
range.  The  epitaxial  layers  were  grown  with  a step-function  change 
in  the  dopant  gas  flow  during  the  continuous  deposition;  during 
the  deposition,  the  dopant  gas  flow  was  changed  from  one  of  the 
above  described  to  the  other  to  simulate  the  step  input.  Both  the 
increasing  and  decreasing  steps  were  used  in  these  experiments. 

The  dopant  profiles  in  the  epitaxial  layers  were  obtained  by 
capacitance-voltage  measurements  on  deep-depletion  MOS  structures 
and  on  planar  and  mesa  p-n  Junctions.  The  thickness  and  dopant 


concentrations  In  the  samples  were  chosen  to  be  compatible  with 
the  C-V  technique.  Spreading-resistance  measurements  and  stacking- 
fault  thickness  determinations  were  used  to  confirm  the  capacitance- 
voltage  measurements.  Figures  13  and  14  are  the  resulting  dopant 
profiles.  To  ensure  that  the  experimental  profiles  were  not  limited 
by  the  resolution  of  the  capacitance-voltage  data-reduction  technique, 
a theoretical  capacitance-voltage  curve  was  generated  by  solving 
Poisson's  equation.  The  data-reduction  program  was  then  applied 
to  this  curve.  The  results  revealed  that  the  experimental  profiles 
varied  slowly  compared  to  the  resolution  profile,  indicating  that 
the  major  features  observed  were  not  artifacts  of  the  analysis 
technique. 
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DISTANCE  FROM  SURFACE 

Fig.  13.  Experimentally  observed  dopant  concentration  as  a function 
of  distance  from  the  surface  of  an  epitaxial  film  for  a decreasing 
step  change  in  the  dopant  gas  flow.  Also  shown  is  an  exponential 
fit  to  the  experimental  curve. 


As  can  be  seen  in  Figures  13  and  14,  the  transition  from  one 
dopant  concentration  to  the  other  occurred  in  approximately  1.2  u» 
corresponding  to  roughly  2 min.  The  heat  cycling  during  the  fabri- 
cation of  the  MOS  and  p-n  Junction  C-V  samples  produced  a square- 
root  of  Dt  of  only  0,07  p because  of  the  low  dlffuslvlty  of  the 
arsenic  dopant  in  the  epitaxial  layer.  Spreading-resistance  data 
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Fig.  14.  Experimental  dopant  profile  for  an  increasing  step  iliange 
in  dopant  gas  flow  and  exponential  fit  to  the  experimental 
prof ile. 


confirm  that  there  is  no  significant  deviation  between  samples 
measured  after  epitaxial  deposition  and  after  complete  C-V  sample 
fabrication.  Tlie  increasing  and  decreasing  steps  reveal  an  abrupt 
onset  of  the  transition  region  followed  by  a gradual  decay.  As 
seen  in  Figures  11  and  14,  the  decay  curve  can  be  fitted  to  an 
exponent)  .1  curve  with  fairly  good  accuracy,  and  a decay  length  of 
0.41  u corresponds  to  a time  constant  of  41  sec.  As  a result, 
the  transient  response  of  the  epitaxial  factor  can  be  approximated 
by  a single  exponential  function,  and  the  response  is  given  by 


H(t)  = 1 


(19) 


where  = 40  sec. 

2.3.3b  SystcMii  analysts.  Altc'r  the  transient  response  of 
a system  has  been  characterized,  it  is  possible  to  predict  the 
output  lor  a certain  input  One  ol  the  simplest  methods  i.s  to 
use  the  convolutional  integral. 


In  a linear  system,  the  output  g(t)  to  an  input  f(t)  can  be 
calculated  by  convolving  the  input  with  the  Impulse  response  l(t). 


(20) 


g(t)  = / 1(t)  f(  t - t)  di 

' r\ 


The  impulse  response  can  be  determined  by  taking  the  time  derivative 
of  the  step  response  h(t);  I(t)  = 9h(t)/9t.  Because  h(t)  = 0 for 
negative  values  of  t,  I(t)  is  also  zero  for  negative  values  of  T. 

The  lower  limit  of  integration  canbe  simplified,  therefore,  to 

t 

g(t)  = ! 1(t)  f(  t - t)  dT  (21) 

0 

The  limits  of  integration  can  be  further  simplified,  depending  on 
the  form  of  the  input.  For  example, if  f(t)  = 0 for  t < 0,  then 
f(t  - t)  = 0 for  t < T,  and  thus  the  upper  limit  of  the  integration 
becomes 


t 

g(t)  = f I(t)  f(t  - t)  dT 
0 


(22) 


Equations  (21)  and  (22)  can  be  evaluated  analytically  or  numerically 
through  numerical  Integration. 

Similar  calculations  can  be  performed  in  the  frequency  domain 
via  the  Laplace  or  Fourier  transform.  The  advantage  of  the  Laplace 
transform  is  that  it  is  possible  to  draw  analogs  to  electrical 
circuits;  the  advantage  of  the  Fourier  transform  is  that  numerical 
methods  are  very  efficient.  In  both  techniques,  repeated  transforms 
are  necessary  to  convert  the  information  from  the  time  domain  to 
the  frequency  domain  and  then  back  to  the  time  domain. 


In  the  epitaxial  reactor,  differentiating  Eq . (19)  obtains 


(23) 


Because  the  growth  rate  G is  constant  in  time,  the  doping  density 
N(x)  is 

N(x)  = G g(t)  (24) 

By  combining  Eqs.  (22),  (23),  and  (24), 

00 

N(x)  • ~ f exp  ( - ^ ) f(t  - 
10  1 ' 


t)  dT 


(25) 


As  noted  earlier,  the  limits  of  integration  can  be  simplified, 
depending  I'li  the  type  of  f(t). 

2.2.3c  Verifioation  of  the  approach.  In  a series  of  exper- 
iments, the  dopant  gas  flow  consisted  of  an  increasing  step  fol- 
lowed by  a decreasing  step,  approximating  a pulse  in  the  input. 

Four  pulse  widths  were  used  (4.3,  3.  1.8,  and  0,8  min).  llie 
resulting  doping  profiles  in  the  substrates  were  measured  using 
the'' C-V  techniques,  and  Figure  If)  plots  the  results.  For  pulse- 
widths  of  4.3  and  3 min,  the  resulting  prf)files  appear  to  be  a 
simple  superposition  of  profiles  obtained  in  Figures  13  and  14 
because  the  pulse  widths  are  longer  compared  to  the  system  time 
constant  (41  sec).  Wlien  using  shorter  pulses  of  1.8  and  0.8  min, 
however,  the  higher  limit  of  doping  was  never  reached,  thereby  in- 
dicating the  limits  of  the  deposition  system  in  responding  to  arbi- 
trary changes  in  the  rate  of  dopant  gas  flow. 

To  verify  the  validity  of  the  transfer-function  approach, 
doping  profiles  were  calculated  using  Eq . (25)  for  the  four  pulse 
inputs.  The  results  are  plotted  in  Figure  15.  There  appears  to 
be  excellent  agreement  between  theory  and  experiments  except  at  the 
onset  of  the  high-to-low  region.  Even  for  the  shortest  pulse, 
where  the  pulse  width  of  48  sec  is  comparable  to  the  system  time 
constant  of  41  sec,  the  model  works  with  good  accuracy. 

In  another  experiment,  the  shape  of  the  dopant  gas-flow  input 
function  was  a ramp  followed  by  a step.  The  results  of  measurements 
and  theoretical  calculations  are  plotted  in  Figure  16  which,  again, 
displays  excellent  agreement  except  for  one  discrepancy  — the 
initial  and  final  values  of  the  dopant  gas  flow  are  the  same  but 
are  different  for  doping  profiles  in  the  epitaxial  layer.  This 
can  be  attributed  to  some  mechanical  malfunction  in  one  of  the 
f 1 owmeters . 


2.4  Summary  for  epitaxy 

This  section  has  defined  several  key  problem  areas  in  modeling 
silicon  epitaxy.  The  kinetics  of  film  growth  and  dopant  Inclusion 
have  been  discussed.  First  order  models  have  been  presented.  In 
the  final  subsection  recent  results  in  transfer  function  modelii.g 
work  in  progress  has  been  discussed. 
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Kig.  16.  Kfspoiise  ot  the  re.iitcr  to  ,i  r;imp  followed  by  a decreasing 
step  input  in  the  dopant  gas  flow. 


3.  SILICON  OXIDATION 
3.1  Introduction 


Modeling  the  kinetics  of  silicon  o.xidation  is  well  known  118]. 

Tlie  first-order  model  is  based  on  the  dll  fusion  and  surface  reaction 
of  a single  oxygen  species  at  the  silicon  intirlace.  Present  trends 
in  the  processing  technology,  iu)wever,  are  quickly  outstripping 
the  utility  of  the  simple  model.  Kor  examp’”,  concentration- 
dependent  effects  can  change  oxide  thickness  by  as  much  as  300% 
for  otherwise  identical  conditions.  l.i>w-temperature  processing 
makes  this  effect  even  more  pronoumed.  Partial  pressures  of 
oxygen,  MCI,  and  other  ambients  also  (hange  growth  parameters. 
Two-dimension.al  local  oxidation  with  inherent  strain-dependent 
effects  is  now  a common  technology  tool  which  requires  more  extensive 
understanding  of  physics  beyond  the  SiO_>-Sl  interface  (41  |.  In 
addition,  there  are  electrical  alterations  ot  the  SIOt-SI  interface 
that  result  from  [)hysical  properties.  Deal  |4Jl  has  recently 
reviewed  this  subject  ;ind  it  is  cle.ir  that  lurther  work  in  this 
area  is  retjuired.  The  presem c oi  oxide  and  surt.u'e  states  can 
substantially  change  MDS  threshold  leakage,  surl.ice  mobility,  and 
overall  device  stability. 

I 

I 

I 


The  mult Istream  models  for  oxidation  kinetics  and  two-dimensional 
dependencies  will  not  be  discussed  here.  It  is  anticipated  that 
model  development  in  these  areas  will  benefit  directly  from  other 
numerical  modeling  work  that  couples  and  conserves  diffusing  particles. 
The  following  discussion  will  be  confined  to  problem  formulations 
that  have  evolved  as  direct  extensions  of  the  first-order  oxidatic’n 
model  [18].  This  model  will  be  reviewed  and  experimental  results 
and  interpretations  will  be  considered  for  orientation,  concentra- 
tion, and  ambient-dependent  oxide  growth  conditions.  In  all  cases, 
modifications  amount  to  coefficient  changes  in  the  basic  model. 


3.2  Tlie  basic  model 

The  first-order  model  solves  for  the  oxide  film  growth  by  integrating 
a flux  equation  for  oxygen  as  it  crosses  the  silicon-oxide  layer  ot 
thickness,  x^,  and  reacts  at  the  silicon  surface.  Figure  17  is  a 
schematic  of  the  model,  with  the  appropriate  flux  terms.  According 
to  Henry's  Law,  the  equilibrium  solid  concentration  is  proportional 
to  the  bulk  gas  partial  pressure,  Pq, 


C = H P 


G 


(26) 


where  C is  the  maximum  silicon-oxide  concentration  for  a given  P(; 
and  H is  the  Henry's  Law  coefficient.  Because  we  are  treating  non- 
equilibrium, the  solid  value  is  less  than  C*.  The  flux,  Fj,  is 
determined  by  the  difference  in  solid  concentrations. 


F^  = h(C  - 


C ) 
o 


(27) 


Here,  C is  the  silicon-oxide  surface  concentration  and  h is  the 
mass-t rHnsfer  coefficient.  The  value  of  h in  the  solid  can  be 
related  to  tlie  gas-transfer  coefficient  by 


h = 


HkT 


(28) 


where  h^  is  the  mass-transfer  coefficient  in  the  gas.  The  second 
and  third  flux  terms  shown  in  Figure  17  represent  diffusive  flux 
(Flck's  Law)  and  the  surface  reaction.  Tlte  diffusive  flux  is  deter- 
mined by  the  concentration  gradient  and  the  effective  diffusivity 

ef  f 


I)  ^ 
eff  Dx 


C - C 
' o 1 


eff\ 


(29) 


F2  = Deff 


F3  = kCi 


F,  = h (C -Cq) 


Ki«.  17.  Boundary  and  Flux  conditions  for  the  Gas-SiO^-Si  system 


In  this  case,  C,  is  tiu'  interface  value  of  the  oxidizing  species 
in  the  siliton  oxide.  The  interface  reaction  rate  is 


where  k is  the  sur 1 ace- reac t ion  rate  constant  for  oxidation. 

Using  the  condition  at  steadv  state  (F|  = F^  = 
solve  tor  C.  anil  C in  terms  ol  C* : 


(31) 


C 


1 


and 


1 


* 

C 


, kx 

h D _ 
erf 


(32) 


Two  limiting  cases  can  be  defined.  Surface-reaction  rate  control 
occurs  when  becomes  large  (compared  to  the  surface  rate  constant 

k)  and 


(33) 


This  condition  implies  that  surface-reaction  kinetics  is  the  slowest, 
and  hence  dominant,  effect.  Mass  transport,  or  diffusion  control, 
occurs^when  k >>  and  , approaches  zero.  In  turn, 

Cq  = Hk  T Cq  where  Cq  is  the  gas-phase  oxidation  concentration. 

Under  these  conditions,  the  flux  toward  the  surface  becomes 


o 


and  the  dominant  effect  is  transport  across  the  oxide  layer. 

To  describe  the  rate  of  oxide  growth,  the  flux  equation  at  the 
interface  can  be  written  as 


dx 

N — ^ = F 
1 dt  3 


kC 


, kx 
^ ®eff 


(35) 


where  Nj  is 
vo 1 ume  (2.2 
ambient) . 
and  t : 


the  number  of  oxidant  molecules  incorporated  per  unit 
X 10^  /cm^  for  dry  O2,  and  4.4  x lO^^/cm^  in  a wet 
The  following  integral  defines  the  relationship  of  x 

o 


+ 


kC*  f dT 
o 


(36) 


wliicli  results  in 


X +Ax  =B(t+i) 
o o 


(37) 


where 


^ ‘^e  f f ( k h ) 


(38) 


B = 


' '^eff  ^ 


(39) 


and 


X . + Ax . 
1 1 


(40) 


The  quadratic  equation  for  x can  be  solved  by  the  following  limit- 
ing cases. 


Case  1 


B A 

x as  - (t  + T)  for  t + T <<  -— 

o A 4B 


(41) 


and  B/A  is  the  "linear"  rate  constant. 
Case  2. 


2 A 

X =B(t+i)  fort+t>>  — 

o 4B 


(42) 


ind  B is  the  "parabolic"  rate  constant. 


Examination  of  the  limiting  forms  indicates  that  the  oxidation 
process  in  the  parabolic  domain  is  diffusion  limited  and,  in  the 
linear  region,  it  is  surface-reaction  limited  (generally,  h >>  k) . 
Both  regions  should  be  directly  dependent  on  the  equilibrium  con- 
centration of  oxidant  in  the  oxide.  Factors  affecting  the  diffu- 
sion process  should  be  most  influential,  therefore,  over  long 
oxidation  times. 


Experiments  have  indicated  that,  with  an  intlally  clean  bare  o 
Si  surface,  Xj^  = 0 for  wet  O2  oxidation,  and  an  effective  200  A 

is  found  for  dry  O2  oxidation  (as  a result  of  an  initial  phase  of 
rapid  oxidation  by  a different  mechanism). 


At  the  present  time  B,  B/A,  and  i are  known  only  for  dis- 
oriented lightly  doped  conditions.  Linder  these  restrict  it)tis , these 
parameters  can  be  expressed  in  the  following  form: 


-E  /kT 

B = C^e 


(4  3) 


-E  /kT 

B/A  = C^e 


(44) 


T = (X^  + AX. )/B 
1 1 


where 

dry  0^: 


wet  O^: 


= 7.72  X 10^  p^/hr 
= 6.23  X 10^  u/hr 

= 28.5  kcal/mole  = 1.23  eV/molecule 
E^  = 46.0  kcal/mole  = 2.0  eV/molecule 

= 2.14  X 10^  p^/hr 
= 8.95  X 10^  p/hr 

Ej^  = 16.3  kcal/mole  = 0.71  eV/molecule 
E^  = 45.3  kcal/mole  = 1.97  eV/molecule 


(45) 


Orientation,  the  addition  of  a chlorine  species,  and  heavy 
substrate  doping  have  all  been  shown  to  affect  B,  B/A,  or  both 
rate ' constants . The  following  sections  discuss  these  individually. 


3.3  Silicon  crystal  orientation  effects 

A log- log  plot  of  oxide  thickness  vs  oxidation  time  for  the  thermal 
oxidation  of  silicon  in  dry  oxygen  at  temperatures  ranging  from 
700°  to  1200°C  is  shown  in  Figure  18.  As  expected,  the  difference 
in  oxidation  rate  between  (100)-  and  (lll)-oriented  silicon  decreases 
gradually  from  700°C  (where  the  differences  are  on  tlie  order  of 
40%)  to  1200°C  (where  the  differences  are  less  than  2%  for  all 
points)  with  the  (111)  oxidizing  faster  than  the  (100).  It  should 
be  noted  that,  in  the  region  below  100  A at  700°C,  the  oxide  thick- 
nesses obtained  for  (100)  and  (111)  silicon  are  essentially  the 
same.  Only  after  the  oxide  thickness  exceeds  100  A do  the  normal 
orientation  differences  appear.  This  observation  suggests  that, 
in  the  initial  oxide  growth  regime  at  700°C,  the  silicon  surface 
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plays  a minimal  role.  Such  conclusions  are  consistent  with  previous 
data  (19 1 obtained  by  Auger  analysis  of  thin  oxides  grown  between 
200°  and  800°C,  which  show  that  the  hydration  of  the  silicon  sur- 
face plays  a crucial  role  in  silicon  oxidation.  The  results  can  be 
of  interest  in  the  fabrication  of  MNOS  devices. 

ilie  (jxidation  rate  data  were  evaluated  using  the  general 
relationship  for  the  thermal  oxidation  of  silicon  given  in  Eq . (37). 
Because  the  parabolic  rate  constant  B depends  on  the  partial  pressure 
of  oxidant  in  the  oxidizing  ambient,  the  solubility  of  the  oxidizing 
species  in  Si02»  and  the  effective  diffusion  coefficient  of  the 
oxidizing  species  in  Si02  (18),  there  should  be  no  variation  of  B 
with  silicon  orientation.  On  the  other  hand,  the  linear  rate  con- 
stant B/A  Involves  the  reaction  rate  constants  at  the  SIO2  surface 
and  at  the  silicon  surface  (18|  so  that  any  orientation  dependence 
of  the  oxidation  rate  should  appear  In  this  term  via  the  constant 
A. 


Because  the  method  of  analysis  utilized  depends  on  the  correction 
factor  T,  which  is  determined  by  extrapolation  and  is  thus  prone 
to  some  error,  oxidation  times  at  least  three  times  greater  than 
the  1 value  for  any  set  of  oxidation  conditions  were  generally  used 
for  the  determination  of  B and  B/A.  in  this  way,  any  error  in  [ 
should  have  little  effect  on  the  calculated  values  of  B and  B/A. 

Tfie  results  of  least-squares  analyses  of  oxidation-rate  data  for 
(100)-  and  ( 1 1 1 )-orienied  silicon  wafers  using  Eq . (37)  are  given 
in  Table  4.  It  should  be  mentioned  that  the  oxide  thickness  (xj) 
at  zero  oxidation  tlmt'  as  determined  by  extrapolation  of  the  x„  vs 
t curve  to  zero  oxidation  time,  was  found  to  be  160  i 40  X for  all 
oxidation  conditions.  Table  4 demonstrates  that,  according  to 
theory,  there  is  essentially  no  difference  in  B for  the  two  orienta- 
tions, at  least  from  1200'’  to  900°C.  Below  900‘’C,  very  long  oxida- 
tion times  (>  100  hr)  become  necessary  to  eliminate  the  effects  of 
:,  and  the  scatter  observed  was  considerable.  At  700°C,  a value 
for  B could  not  be  obtained  because  a plot  of  oxide  thickness  versus 
time  yielded  a straight  line,  indicating  that  the  parabolic  contri- 
bution to  the  overall  oxidation  rate  was  negligible.  The  linear 
rate  constant  B/A,  shows  the  expected  orientation  dependence,  with 
B/A  for  (111)  always  being  equal  to  or  greater  than  that  for  (100). 

It  should  be  noted  that  the  kinetic  parameters  for  (11 l)-oriented 
silicon  in  Table  4 are  in  very  close  agreement  with  those  obtained 
by  Deal  and  Grove  [18].  Such  agreement  is  remarkable  when  one 
considers  the  different  oxidation  conditions  (humidity,  gas  quality, 
flow  rates,  furnaces,  etc.)  prevailing  in  iwo  laboratories  carrying 
ut  similar  experiments  11  years  apart. 

Activation  energies  were  obtained  by  fitting  the  data  in  Table 
4 to  an  Arrhenius  equation  in  the  form  of 


k = k e 
o 


-E  /kt 
a 


where  Eg  is  the  activation  energy,  T is  the  temperature  in  ‘’K,  and 
k is  a gas  constant.  Plots  of  B and  B/A  vs  temperature  are  shown 
in  Figures  31  and  32  in  connection  with  the  HCl  data  described 
later.  This  analysis  obtained  values  of  1.3  eV  (30  kcal/mole)  and 
1.2  eV  (28  kcal/mole)  for  the  activation  energy  of  B for  (100)- 
and  (lll)-oriented  silicon,  respectively,  which  agrees  well  with 
previously  determined  values  for  the  dry  thermal  oxidation  of 
( 1 ll)-orlented  silicon  [19,  20]. 

A similar  analysis  for  the  linear  rate  constant  for  the  (111) 
orientation  yielded  a value  of  2 eV  (47  kcal/mole),  which  is  in 
close  agreement  with  the  energy  required  to  break  a Si-Si  bond 
[21].  Determination  of  the  activation  energy  for  B/A  for  (100)- 
oriented  silicon,  however,  affords  two  choices.  if  the  data  from 
900°  to  1200°C  are  analyzed  and  the  low-temperature  points  are 
neglected  because  of  large  scatter,  a value  of  2.5  eV  (57  kcal/mole) 
is  obtained  — an  increase  of  ~ 25%  over  that  for  (100).  This 


RATE  CONSTANTS  FOR  <111>  AND  -^100>  SILICON 


u 

x: 


(N  GO 

\ LO  r- < ir>  fN  rr 

E r^r^r^oor^ror'jm 

3 voa>vx5Chr^r^f-Hmoo 

TT  00  TT  r-'  o»— I o o oo 

< »HiH  <N(N  OO  OO  OO  OO 

\ 

CQ  r-H  rH  OO  OO  OO  OO  OO 


u ir> 

x:  m VO  CM  rH  in 

\ moo  r^Tf  iHin  moo  r^rH 

04  Lin  oo  mm  oh 

g (NCM  OO  OO  OO 

a oo  oo  oo  oo  oo 

II 

cn  oo  oo  oo  oo  oo  i i 


m 

^ ov  O 

•H 

m o 

»H 

E m Tj- 

O 00 

av  <N 

CM  iH 

m 

;i  O O 

iH  O 

H ^ 

^ <N 

m 

• • 

• t 

• • 

• • 

• 

• 

< o o 

O O 

o o 

o o 

o 

o 

u mm  <T»  o>  mm 

^ oo  OO  mm  (n<n  oo 

• • • • • • • • • • I I 

H oo  oo  oo  mH  om  it 


Increase  was  predicted  by  Ligenza  from  consideration  ot  steric 
hindrance  in  the  (100)  and  (111)  planes  of  silicon  during  thermal 
oxidation  [ 22  ] . 

On  the  other  hand.  If  the  data  from  700°  to  1200°C  are  analyzed, 
the  900°,  1100°,  and  1200°C  points  could  be  considered  as  scatter 
from  a line  drawn  parallel  to  the  one  for  (111)  orientation,  passing 
through  the  700°,  800°,  and  1000°C  points,  thereby  resulting  in 
the  same  activation  energy  as  (111 )-or iented  silicon.  Indeed, 
some  faith  could  be  put  in  such  as  assumption,  because  the  700"C 
data  give  B/A  directly  from  an  vs  t plot  as  indicated  above. 

One  could  argue,  therfore,  that  the  higher  temperature  analyses  were 
in  error  due  to  the  inaccuracy  in  determining  A from  a plot  of  x^, 
vs  (t  + t/Xj^).  The  slopes  of  these  lines  (B)  are  quite  large; 
therefor',  a small  change  in  the  slope  would  shift  the  Intercept 
(A),  and  thus  B/A,  considerably.  At  present,  it  is  not  known  which 
of  the  above  alternatives  is  the  correct  one.  In  any  case,  the 
preceding  discussion  indicates  that  there  has  been  little  or  no 
work  reported  since  19A1  concerning  the  differences  between  oxid.a- 
tion  rates  and  corresponding  activation  energies  of  (100)-  and 
( 1 1 1 )-oriented  silicon  in  dry  oxygen. 


3. A Impurity  doping  effects 

Tile  effects  of  Impurity  doping  levels  on  thermal  oxidation  rates 
are  intimately  connected  with  a widely  encountered  phenomenon  in 
semiconductor  processing  — namely,  Impurity  redistribution. 

Figure  19  Illustrates  how  redistribution  and  thermal  oxidation 
interact . 

As  a thermal  oxide  is  grown  over  a doped  silicon  substrate, 
redistribution  of  the  impurity  results.  In  the  case  of  phosphorus, 
arsenic,  and  antimony,  the  dopant  atoms  tend  to  pile  up  at  the 
surface,  resulting  in  a higher  surface  concentration  than  background 
concentration  (C^  >>  Cg) . In  boron,  the  opposite  effect  takes 
place,  resulting  in  surface  depletion  (Cg  « Cg) . 

Very  heavily  doped  substrates  (Cg  >-  10^*^),  it  has  been  observed 
[23]  that  for  both  phosphorus  and  boron,  the  oxidation  rates  can 
be  substantially  different  (generally  faster)  than  those  observed 
on  lightly  or  moderately  doped  substrates.  With  respect  to  Figure 
19,  the  two  parameters  that  have  been  correlated  [23]  with  this 
increased  oxidation  are  Cg  (the  dopant  concentration  in  the  silicon 
at  the  surface)  and  C^j^^  (the  average  impurity  concentration  in  the 
oxide ) . 


1 


Fig.  19.  Effect  of  redistribution  on  impurity  doping  profiles. 


Intuitively,  tlie  effect  of  is  to  reduce  the  amount  of  energy 
required  to  break  .Si  bonds  and  thus  to  affect  the  surface  reaction; 
it  would  thus  be  expected  to  influence  the  linear  rate  constant  B/A 
and  not  B.  The  effect  of  is  to  change  the  diffusiozi  constant 

for  the  (>2  or  11)0  oxidizing  species  in  the  Si02,  nnd  it  would  thus 
be  expected  to  affect  the  parabolic  rate  constant  B. 

To  study  these  enhanced  oxidation  rates  in  more  detail,  a 
large  numbc’r  of  heavily  doped  wafers  were  prepared  as  illustrated 
in  Figure  20.  Because  of  the  difficulty  in  purchasing  wafers  with 
doping  levi’ls  ■ 10*'^/cm^,  the  samples  were  prepared  by  diffusion 
ol  phosphorus  into  standard  lightly  doped  (111)  substrates.  The 
If  predeposition  and  drive-in  schedules  were  chosen  so  that  the  diftused 

proliles  were  flat  (within  10%)  over  the  first  2 p into  tlie  silicon. 
I'he  drive-ins  were  done  in  N2  to  prevent  oxidation  and  hence  redistri- 
bution during  the  form.ition  of  heavily  doped  layers. 

Five  types  of  s.amples  were  prepared  with  surface  doping  levels 
between  9 x 10*'^/cm^  and  solid  solubility  S x lO^^Vcm^),  as 
indicatt'd  in  Figure  20.  Tlie  resulting  diffused  profiles  were 
measured  using  spreading  resistance  and  anodic  sectioning  techniques. 
Agri'ement  between  the  two  measurement  techniques  was  excellent. 


It  is  important  to  note  that  the  doping  concentrations  tabulated 
in  Figure  20  represent  electrically  active  concentrations,  not 
chemical  concentrations. 
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Fig.  20.  Surface  concentrations  of  diffused  heavily  doped  samples. 


Following  preparation  of  the  samples,  a series  of  initial 
experiments  were  conducted  to  evaluate  the  itwigii i tude  of  the  enhanced 
oxidation  effect  over  regions.  A typical  result  is  shown  in 
Figure  21.  The  horizontal  scale  is  again  an  electrically  active 
surface  concentration.  The  data  indicate  appro xinvitely  a three- 
to-one  range  in  oxide  thickness  for  these  wafers  which  were 
oxidized  simultaneously. 

The  choice  of  a low  oxidizing  temperature  tends  to  mtiximize 
the  enhanced  oxidation  rate  for  the  common  N-tvpe  Impurities  because 
as  explained  above,  they  tend  to  pile  up  at  the  silicon  surface. 

As  a result,  the  dominant  effect  is  on  the  reaction  occurring  at 
the  S1-S102  Interface,  and  hence  the  linear  rate  constant  B/A 
which  dominates  the  overall  reaction  at  low  temperatures.  It  should 
be  noted  that  low-temperature  oxidation  is  becoming  Increasingly 
important  in  the  semiconductor  industry  because  of  the  trend 
toward  larger  diameter  silicon  wafers. 


Fifi.  21.  Oxide  thickness  vs  doping.  Samples  oxidized 
s imul taneous ly  . 


Wet  oxidation  also  tends  to  increase  the  effect  observed  iii 
Figure  21  because  there  is  more  redistribution  for  wet  oxides 
when  N-type  impurities  are  involved  (Cs/Cb  is  larger  for  wet  tlian 
dry  oxides)  |23,24]. 

Following  these  Initial  experiments,  a more  carefully  designed 
series  of  oxidations  were  made  to  evaluate  tlie  enhanced  oxida- 
tion over  N+  regions  under  a wide  variety  of  temperatures  and 
ambient  conditions.  Typical  of  the  results  obtained  are  the  data 
shown  in  Figures  22  and  23.  The  bottom  curve  L\)  in  each  case  is 
lightly  doped  (1  x lO^^)  material,  and  the  results  agree  witli 
previously  published  data  [18].  Tlie  upper  curve  (F)  corresponds 
to  a substrate  doped  to  approximately  solid  solubility. 

The  following  observations  can  be  made  about  tliase  data,  even 
in  the  form  shown  in  these  figures. 

(a)  The  enhanced  oxidation  effect  is  more  pronounced  at 
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lower  temiH'rat  iires  (900  C)  than  at  higher  temperat  urts 
( 1000°C) . 

(2)  riu'se  effeits  are  also  more  pronounced  for  shorter  timis 
and  thinner  oxides  than  they  are  for  longer  times  and 
tliicker  oxides. 

Both  of  these  observations  are  in  agreement  with  our  previously 
stated  expectations;  that  is,  N-type  dopants  that  pile  up  at  the 
silicon  surface  and  segregate  into  the  Si  should  affect  the  reaction 
kinetics  at  the  Si-Si02  interface  far  more  than  they  affect  tlie 
diffusion  of  the  oxidizing  species  through  the  Si02.  As  a result, 
we  would  expect  B/A  to  be  affected  and  not  B,  in  agreement  witli 
the  observations  stated  above. 

Based  on  the  general  oxidation  relationship,  we  can  consider 
from  a theoretical  point  of  view  what  the  curves  in  Figures  22  and 
2 3 should  look  like  if  the  only  effect  of  tr*"  doping  was  on  the 
surface-reaction  kinetics.  Figure  24  is  a series  of  theoretical 
curves  generated  from  this  oxidation  relationship.  The  various 
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Fig.  24.  Theoretical  curves  of  oxide  thickness  vs  time  for 
various  linear  rate  constants. 
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curves  were  generated  for  values  of  the  linear  rate  constant  B/A 
between  the  normal  or  lightly  doped  substrate  value  (2.8  x 10 
U/mtn  at  900°C)  and  500  times  this  value  (corresponding  to  a much 
enhanced  surface-reaction  rate).  The  qualitative  agreement  between 
these  theoretical  curves  and  the  experimental  results  in  Figure 
22  is  apparent . 


The  analyses  of  these  data  in  Figures  22  and  23  can  be  carried 
one  step  further  by  extracting  linear  and  parabolic  rate  constants 
according  to  standard  techniques  [18].  The  results  for  the  900°C 
oxidation  are  shown  in  Figure  25.  The  following  conclusions  can 
be  drawn  from  these  data. 

(1)  The  parabolic,  rate  constant  is  essentially  unchanged 
by  the  N"*"  substrate.  This  is  in  accordance  with  our 
expectations  because  most  of  the  phosphorus  segregates 
into  the  silicon. 

(2)  Tlie  linear  rate  constant  is  virtually  unchanged  for 
doping  levels  ^ 5 x 10^^/cm^.  Above  this  concentration 
B/A  increases  approximately  one  order  of  magnitude  as 
the  doping  level  is  increased  to  the  so  1 id-sol ub i 1 i tv 

1 imit . 


lO'^  10'*  I0‘*  10*®  • 10*' 

PHOSPHORUS  CHEMICAL  CONCENTRATION 


Fig.  25.  Linear  and  parabolic  rate  constants  vs  chemical  dopant 
concentrations  at  900 *C. 
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It  slimild  hi-  nuti'd  that  the  horizont.il  scale  in  this  fiKure 
h.is  been  I'orn-i’ted  to  .i  chem ic.i  1 phospliorus  concentration  in  the 
siihstrate,  r.ither  than  the  e 1 ec  t r i la  1 1 y active  concentration 
previously  used.  The  chemical  phosphorus  concent r. it  ion  (and  hence 
the  pi-rcenta>;e  of  tlie  species  electrically  active  in  our  d i f 1 used 
s.iniples)  was  determined  by  sputtering  Au>;er  teclini<|ues  by  Joh.mnesen 
et  .il.  .It  St.inford  University  [25  1.  There  is  still  some  question 
icMicernin^  the  absolute  accuracy  ot  the  horizontal  scale  in  FiKure 
J)  bee  .luse  of  the  sensitivity  limitations  of  the  Au>»er  technique. 
Independent  me.isurements  on  these  samples  are  currently  being 
performi-d  using  SIMS  techniques.  These  should  resolve  the  question 
ol  chemii'.il  vs  electrically  .active  phosphorus  concen  t r.a  t i ons  in 
these  s.iniples. 

The  10()()°U  e.'ctracted  r.ite  constants  show  a form  very  simil.ir 
to  Figure  25.  Samples  have  also  been  oxidized  at  1100°  and  800°C, 
•ind  the  resulting  oxide  thicknesses  are  being  measured.  Completion 
of  thi-se  d.it.i  at  the  four  temperatures  will  allow  us  to  plot  both  B 
.ind  B/A  .IS  .1  function  of  temperature  .and  to  extract  activation 
energies  tor  the  rate  constants  .as  a function  of  substrate  doping 
1 evel  . 


The  re|)orted  activ.ation  energy  for  lightly  doped  w.afer  for  B/A, 
the  linear  rate  constant,  is  approximate  1 v 46  kcal/mole  1181. 

This  comp.ires  well  to  the  energy  required  to  break  a Si-Si  bond 
which  is  42.2  kcal/mole  [21  ].  The  900°  and  1000°C  data  indic.ate 
tii.it  B/A  becomes  extremely  large  for  the  samples  doped  close  to 
solid  solubility.  The  activation  energy  for  B/A  in  this  case, 
tiiereiore,  .ipproailies  zero,  thereby  indicating  that  tlie  re.action 
.It  tile  Si-SiOj  interf.ace  is  no  longer  a rate-limiting  step  at 
these  doping  levels;  that  is,  it  appears  that  the  incorporation  of 
di  |).int  .itoms  in  concentrations  close  to  solid  solubility  so  strains 
the  l.ittice  structure  tliat  bond  breaking  becomes  extremely  easy. 
Wlietiier  this  simplified  model  is  correct  can  only  be  verified  by 
.1  better  model  ot  tile  molecular  and  atomic  structure  near  the  Si- 
Si()2  interlace.  Sucli  a model  is  necessary  also  for  prediction  of 
c h.irge  densities  and  is,  therefore,  a goal  of  our  research  progr.im. 
We  .in-  .also  beginning  to  pursue  it  in  the  manner  indicated  at  tlie 
beginning  of  this  section. 


1.5  Oxid.it  ion  in  HCl/O^  mixtures 

One  ol  the  most  significant  developments  in  the  passivation  of 
thermally  grown  silicon-dioxide  films  over  the  past  several  years 
has  been  the  addition  of  a chlorine  species  during  silicon  oxida- 
tion. Il  has  been  demonstrated  (26,27,28]  that  such  additions 

result  in  improved  threshold  stability  and  increased  dielectric 
strength.  In  addition,  it  has  been  observed  that  chlorine  additions 
increase  the  rate  of  silicon  oxidation  [29,30,31],  but  extensive 
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data  for  widely  varying  oxidation  conditions  have  not  been  reported. 

In  this  investigation,  therefore,  the  thermal-oxidation  kinetics 
of  silicon  in  O2/HCI  mixtures  is  cliaracterlzed  as  a function  of 
oxidation  temperature,  HCl  concentration,  and  silicon  orientation; 
the  results  are  compared  to  those  obtained  from  silicon  oxidations 
performed  in  a dry-oxygen  atmosphere. 

3.5.1  Rate  constants.  Log-log  plots  of  oxide  thickness  vs 
oxidation  time  for  the  thermal  oxidation  of  (100)-  and  (111)- 
oriented  silicon  in  O2/HCI  mixtures  at  temperatures  of  900’, 

1000°,  and  1100°C  are  shown  in  Figures  26,  27,  and  28,  respectively. 

It  is  clean  that  the  (111)  silicon  always  oxidizes  faster  than  does 
the  (100).  It  can  also  be  seen  that  the  effect  of  the  HCl  addition 
is  to  increase  the  oxidation  rate  relative  to  oxidation  in  dry 
oxygen.  In  particular,  a relatively  large  Increase  in  the  oxidation 
rate  occurs  with  the  initial  1%  HCl  addition.  Subsequent  additions 
have  a somewhat  smaller  effect,  but  the  overall  rate  is  systematically 
increased.  It  can  also  be  observed  in  these  figures  that  the  over- 
all increase  in  oxide  thickness  for  a particular  oxidation  time  on 
increasing  the  HCl  concentration  from  0 to  10%  is  larger  for  (100)- 
than  for  ( 1 1 1 ) -or  iented  silicon.  This  effect  may  be  related  to  the 
silicon-etching  phenomenon.  It  should  be  noted  that  very  few  data 


Fig.  26.  Oxide  thickness 
vs  oxidation  time  for 
the  oxidation  of  <100>- 
and  <lll>-erlented  N- 
type  silicon  in  various 
OgHCl  mixtures  at  900 *C 
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Kig.  27.  Oxide  thickness  vs  oxidation  time  for  the  oxidation  of 
•:lil>  and  <100>  oriented  N-type  silicon  in  various  O2/HCI 
mixtures  at  1000'’C. 


Fig  28.  Oxide  thickness 
vs  oxidation  time  lor 
the  oxidation  of  111  - 
and  loo  -oriented  N- 
type  silicon  In  various 
O2HCI  mixtures  at  1100 *C. 
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points  are  plotted  for  oxidations  carried  out  at  1100°C  with  HCl 
concentrations  greater  than  1%  because  of  silicon  etciilng  and 
oxide-bubble  forming  oxidation  with  O2/HCI  mixtures  at  1100°C 
[30,31,32,33] . 

To  separate  the  effects  of  HCl  on  the  parabolic  B and  linear 
B/A  rate  constants  and  to  obtain  some  indication  of  the  role  HCL 
plays  1'  silicon  oxidation,  Eq.  (37)  was  utilized.  During  the 
determination  of  i,  Xj^  was  assumed  to  be  essentially  constant  at 
160  ± AO  A for  all  oxidation  conditions  investigated.  Although 
xj^  should  decrease  with  HCl  addition  due  to  water-vapor  generation 
as  per  the  reaction  of  O2  and  HCl 

A HCl  + 0^  -*■  2 H^O  + 2 Cl^  (‘^7) 

the  extrapolations  performed  to  determine  X|  and  thus  x could  not 
detect  small  changes  in  Xj^ . Nevertheless,  because  1 depends  on  Xj, 

B,  and  B/A,  a decrease  in  X was  observed  on  HCl  addition,  probably 
in  part  as  a result  of  water  generation  because  Xj^  and,  therefore, 

I are  zero  for  oxidation  in  a steam  amblent[18]. 

The  results  of  least-squares  analyses  on  the  data  in  Figures 
26,  27,  and  28  are  summarized  in  Figures  29  and  30  . which  are  semi  log 
plots  of  the  effect  of  HCl  concentration  on  the  parabolic  and  linear 
rate  constants,  respectively.  It  should  be  noted  that  the  rate 
constants  generated  via  Eq . (37)  and  plotted  in  Figures  29  and  30  are 
ef^ec_tj^y£  rate  constants  in  that  they  represent  the  combined  effects 
of  oxygen,  water,  and  chlorine  species  on  B and  B/A. 

Figure  29  shows  that  essentially  no  orientation  effect  is 
observed  In  the  parabolic  rate  constant  B,  in  accordance  with  theory; 
also,  increasing  the  HCl  concentration  above  1%  results  in  a linear 
increase  in  B (the  plots  are  put  on  a semilog  scale  for  convenience). 
The  increase  in  B going  from  0 to  1%  HCl  in  dry  oxygen  is  not  as 
consistent  among  the  three  temperatures,  however.  For  1000°  and 
1100°C,  a large  increase  in  B due  to  a 1%  HCl  addition  occurs, 
whereas  this  trend  is  not  observed  at  900°C.  Tliese  observations  are 
in  agreement  with  previous  work  that  revealed  a linear  Increase  of 
B with  a 0 to  9%  HCl  addition  at  900°C  [31]  and  a large  increase  in 
B due  to  a 1%  HCl  addition  at  1100°C  [30].  Presumably,  the  large 
Increase  in  B and  the  linear  Increase  with  subsequent  HCl  additions 
are  the  result  of  water  generation  as  described  in  Eq.  (A7)  because 
small  amounts  of  water  Increase  the  thermal-oxidation  rate  of  silicon 
in  dry  oxygen  [3A],  As  pointed  out  by  van  der  Meulen  and  Cahill  [31], 
however,  gas-phase  thermodynamic  equilibrium  calculations  Indicate 
that  the  amount  of  water  generated  via  Eq.  (A7)  cannot  account  for 
the  observed  Increase  in  the  oxidation  rate.  This  is  consistent 
with  the  published  results  that  a CI2  addition  to  a dry-oxygen 
atmosphere  yields  a considerably  higher  slllcon-oxldatlon  rate  than 


Kig.  29.  Parabolic  rate  constant  vs  % HCl  for  <111>  and  ^ 100> 
oriented  N-type  silicon  at  900°,  1000°,  and  1100°C. 


Fig.  30.  Linear  rate 

constant  vs  percent  HCl 
for  <111>-  and  <100>- 
oriented  N-type  silicon 
at  900*,  1000*,  and  llOO'C. 


an  equivalent  HCl  caddltlon  at  1150°C  [34],  suggesting  tliat  the  chlorine 
species,  rather  than  the  water  generated,  may  be  primarily  responsible 
for  the  increase  in  the  oxidation  rate  at  high  temperatures.  Because 
it  is  difficult  to  obtain  CI2  which  contains  less  than  3 ppm  water, 
however,  some  effect,  although  small,  would  be  anticipated  from  the 
water  present. 

In  regard  to  the  above  argument  concerning  the  amount  of  water 
generated  via  Eq.  (47),  thermodynamic  calculations  such  as  those  in 
Ref.  [31]  give  only  the  Idealized  gas-phase  equilibrium  conditions. 

The  presence  of  a surface  — especially  one  such  as  Si02  which  is 
useful  for  catalytic  purposes  — could  shift  the  equilibrium  and  thus 
alter  the  amount  of  water  generated  and  Incorporated  into  the  oxide 
surface. 

Some  insight  can  be  gained  into  the  silicon-oxidation  regime 
affected  by  water  or  chlorine  [as  generated  via  Eq.  (47)]  from 
consideration  of  the  definitions  of  the  parabolic  (B)  and  the  linear 
(B/A)  rate  constants.  The  definitions  are  as  follows  [18]: 
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where  D ^ is  the  effective 
concentration  of  oxidant  in 
molecules  incorporated  into 
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diffusion  coefficient,  C is  the  imiximum 
oxide,  and  is  the  number  of  oxidant 
a unit  volume  of  oxide. 


(49) 


where  k.  is  the  surface-reaction  rate  constant,  and  h is  the  gas-phase 
transport  coefficient. 


Inspection  of  these  definitions  and  a comparison  of  Figures 
29  and  30  reveal  that  the  monatonic  increase  in  B with  HCl  concen- 
tration is  related  to  the  effective  diffusion  coefficient  Dgff 
because  this  term  does  not  appear  in  the  expression  for  B/A.  It 
also  appears  that  the  rapid  initial  increase  in  B and  B/A  may  be 
related  to  the  solubility  of  the  oxidizing  species  in  the  oxide  film. 
Although  the  latter  conclusion  could  be  explained  by  water  addition, 
this  explanation  would  not  account  for  the  Increase  in  because 

this  parameter  is  known  to  decrease  with  water  addition  [18]. 


In  light  of  the  above  arguments,  the  gradual  increase  of  the 
parabolic  rate  constant  with  increasing  HCl  concentration  could 
be  partially  related  to  a chlorine  species  because  it  has  been 


osl.ib  1 islied  Lliat  a higher  concentration  of  chlorine  is  incorporated 
into  the  silicon-dioxide  film  as  the  HCl  concentration  and/or  the 
silicon-oxidation  temperature  are  increased  [35].  Consequently, 
increasing  the  chlorine  concentration  may  cause  the  Si()2  lattice  to 
be  strained  (especially  near  the  interface)  thereby  allowing  diffusion 
of  the  oxidant  to  occur  more  easily  and  increasing  the  oxidation 
rate  as  observed  in  the  oxidation  of  heavily  boron-doped  silicon  [ 36 j . 
Because  the  ionic  radius  of  boron  when  substituted  for  silicon  in  a 
silicon  lattice  (which  should  be  close  to^the  radius  of  boron 
substituted  for  silicon  in  SIO2)  is  0.88  A [37),  and  the  ionic 
radius  of  chlorine  is  1.8  A [38],  it  would  be  expected  that  the  SiO^ 
lattice  expansion  would  be  even  greater  for  chlorine  as  compared  to 
boron.  This  speculation  is  also  consistent  with  the  previous 
observation  that  CI2  additions  result  in  larger  silicon-oxidation 
rates  than  do  HCl  additions  at  high  temperatures  because  more  chlorine 
is  apparently  incorporated  into  the  SiO->  with  O7/CI2  than  with 
O2/HCI  mixtures  [35]. 

figure  30  reveals  a strong  orientation  effect  on  the  linear 
r.ite  const.ant  B/A,  which  apparently  decreases  with  increasing 
temperature.  Tills  observation  was  discussed  in  the  previous  section 
on  orientation  effects.  As  generally  found  for  B,  an  initial  rapid 
increase  in  B/A  can  also  be  observed  when  increasing  the  HCl  concen- 
tration in  the  oxygen  atmosphere  from  0 to  1%;  however,  a further 
increase  resulted  in  essentially  no  change  in  B/A.  Although  an 
explanation  for  these  observations  is  not  apparent,  it  is  conceivable 
that  the  breaking  of  Si-Si  bonds  could  be  promoted  by  etching  the 
silicon  substrate  which  is  known  to  occur  with  the  use  of  HCl 
[30,31,32,33]  to  increase  the  linear  rate  constant.  Increased 
HCl  concentrations,  however,  should  then  further  increase  B/A  unless 
an  equilibrium  between  oxidation  and  etching  is  established. 

Additional  investigation  of  this  phenomenon  is  required  if  we 
are  to  understand  the  mechanism  Involved.  Preliminary  results 
obtained  by  Deal  with  O2/CI2  oxidation  atmospheres  in  the  Fairchild 
Laboratory  substantiate  the  above  interpretation;  relatively  large 
increases  in  B/A  c.iused  by  a CI2  addition  at  1100°C  were  observed, 
with  tons  iderab 1 y smaller  increases  noted  for  B. 

3.5.2  Activation  energies.  Figure  31  is  an  Arrhenius  plot  of 
the  p.trabolic  rate  constant  B for  the  thermal  oxidation  of  silicon 
in  (>2  and  several  O2/HCI  mixtures.  Because  essentially  no  orienta- 
tion effect  was  observed  for  B,  only  the  (111)  orientation  is 
plot  ted . 


A leas  I -squares  analysis  of  the  data  for  dry  oxygen  produced 
a value  of  1.2  eV  (28  kcal/mole),  which  agrees  well  with  previously 
determined  values  for  the  thermal  oxidation  of  silicon  in  dry 
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Fig,,  il . Arrhenius  plots  of  the  parabolic  rate  constants  for 
silicon  oxidized  in  various  O^/HCl  mixtures.  Silicon 
orientation  is  <111>. 


oxygen  [18,  39]  and  for  the  diffusivity  of  oxygen  through  fixed 
silica  [40].  It  can  be  seen  that  the  HCl  addition  to  the  dry- 
oxygen  atmospheres  resulted  in  a convex  curvature  (relative  to  dry- 
oxygen  oxidation)  of  the  Arrhenius  plot. 

Figure  32  is  an  Arrhenius  plot  of  the  linear  rate  constant  B/A 
for  the  thermal  oxidation  of  silicon  in  O2  and  several  O2/HCI 
mixtures.  A least-squares  analysis  of  the  data  for  dry-oxygen  (0%  HCl) 
oxidation  of  (lll)-orlented  silicon  yielded  a value  of  2 eV  (47  kcal/ 
mole),  which  is  in  close  agreement  with  the  energy  required  to  break 
a Si-Si  bond  [21].  Determination  of  the  activation  energy  for 
(lOO)-orlented  silicon,  however,  affords  two  choices  as  explained 
in  the  previous  section  on  orientation  effects. 


Similar  to  the  parabolic  rate  constant,  the  Arrhenius  plot  of 
the  linear  rate  constant  for  HCI/O2  mixtures  (Figure  32)  is  curved, 
and  the  curvature  is  concave  relative  to  dry-oxygen  oxidation. 

Baaed  on  Eq . (47),  it  is  not  surprising  that  Arrhenius  plots  of  the 
rate  constants  for  oxides  grown  in  HCl  atmospheres  should  result  in 
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Fig.  52.  ArrlicMiius  plots  of  tlie  linear  rate  constants  for 
and  ■'100'*  oriented  silicon  oxidized  in  various  O^/llCt 


■ I 1 1 ■ 
mixtures . 


curved  lines.  it  is  apparent  tiiat  more  than  one  activ.ilion  process 
is  occurring  (such  as  oxidation  caused  by  oxygiui  and  w.iter  vapor 
and  several  possible  effects  resulting  from  chlorine).  As  the 
oxidatictn  temperature  is  changed,  the  relative  import.. ..ce  of  tliesi' 
activ.ated  processes  may  also  cliange.  It  becomi'S  neiess.iry,  I tic  ri' f o r<’ , 
to  separate  out  the  effects  of  water  and  chlorine  on  tlu'  silii-on 
oxidat’on  rate.  Toward  this  end,  experiments  are  being  pt’rformed 
-o  to  investigate  the  oxidation  rate  of  silicon  in  0.,/H,()  mixtures. 

4.  CO  MCI.  LIS  ION 

In  the  previous  section  the  fundamentals  of  silicon  oxidation 
and  first  order  process  models  have  been  presented.  The  subsequent 
discussions  Illustrate  changes  In  the  model  coefficients  which  can 
account  for  changes  In  oxidation  rates  based  on  orientation, 
ambient  and  surface  doping  effects. 


In  both  the  epitaxy  and  oxidation  sections  the  discussion 
has  focused  on  two  areas.  First,  the  basic  process  models  have 
been  defined.  Second,  recent  advances  have  been  summarized  which 
extend  the  range  of  application  of  the  basic  models  to  Include 
state-of-the-art  processing  techniques.  This  includes  programmable 
doping  of  epitaxial  layers  and  ambient  and  surface  effects  during 
oxidation. 
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